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Abstract 



Existence and uniqueness are investigated for a nonlinear diffusion problem of 
phase-field type, consisting of a parabolic system of two partial differential equa- 
tions, complemented by Neumann homogeneous boundary conditions and initial 
conditions. This system aims to model two-species phase segregation on an atomic 
lattice [19]; in the balance equations of microforces and microenergy, the two un- 
knowns are the order parameter p and the chemical potential p. A simpler version 
of the same system has recently been discussed in [8j. In this paper, a fairly more 
general phase-field equation for p is coupled with a genuinely nonlinear diffusion 
equation for p. The existence of a global- in-time solution is proved with the help of 
suitable a priori estimates. In the case of a constant atom mobility, a new and rather 
unusual uniqueness proof is given, based on a suitable combination of variables. 
Key words: phase-field model, nonlinear laws, existence of solutions, new unique- 
ness proof 
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1 Introduction 



In this paper, the last so far in a series [TJ, [HJ [9j Ell E] , we further our mathematical 
analysis of a mechanical model proposed by one of us [19] for phase segregation through 
atom rearrangement on a lattice. On postponing a detailed presentation of the model and 
its antecedents until next section, we begin by pointing out what features of the system 
we study are more general, and therefore more difficult to handle mathematically, than 
in our previous paper [8]. 

The initial and boundary value problem we here tackle consists in looking for two 
fields, the chemical potential \i > and the order parameter p £ (0, 1), solving 

2h(p)dtfi + fih'(p)d t p-dxv(K(fi)V(i) =0 inftx(0,T), (1.1) 

5d tP -Ap + f'(p) = ph'(p), in ft x (0, T), (1.2) 

(k(/i)V/i) • u\ r = d u p\ r = on T x (0, T), (1.3) 

/i(-,0) = /i and p(-,0)=p in ft, (1.4) 

where ft denotes a bounded domain of M 3 with conveniently smooth boundary T, T > 0, 
and d v denotes differentiation in the direction of the outward normal v. In fll.ip . the 
atom mobility is specified by a nonnegative, continuous and bounded, nonlinear function 
k of p (in particular, the degeneracy of k around the critical value p = is admitted). 
The problem is parameterized by two scalar-valued functions, h and /, and two positive 
numbers, e and 5, both intended to be small. The parameter functions enter into the 
definition of the system's free energy 

V = ftp, Vp,p) = -ph(p) + /(p) + ||V P | 2 , (1.5) 

where the last two terms favor phase segregation, the former because it introduces local 
energy minima and the latter because it penalizes spatial changes of the order parameter 
(we have set equal to 1 the relative material constant). For h, one can take any smooth 
function provided it is bounded from below by a positive constant: 

Hp) > |; (1.6) 

for /, the sum 

f{p)=h{p) + h{p) 

of a convex and lower semicontinuous function / x , with proper domain D(fi) C R, and of 
a smooth function f 2 with no convexity properties, so as to allow for a double or multi-well 
potential /. Note that f\ need not be differentiable in its domain, so that its possibly 
multivalued subdifferential (3 := df\ may appear in (11.21) in place of /{; in general, (3 is 
only a graph, not necessarily a function, and it may include vertical (and horizontal) lines, 
as for example when 

f if < p < 1 
h{p ) = l m (p = \ - (1-7) 

[ +oo elsewhere 

and (3 = <9/[ 0i i] is specified by 

£ £ j3(p) if and only if £ 
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The simpler situation dealt with in [5] obtains for k constant-valued (and hence set 
equal to 1, without any loss of generality), 

Hp) = p, (1.9) 

and / a double-well potential defined in (0,1), whose derivative /' is singular at the 
endpoints p = and p = 1: e.g., 

f(p) = a, {p hx{p) + (1 - p) ln(l - p)} + a 2 p (1 - p) (1.10) 

for some positive constants ai and a 2 Q Under these less general circumstances, system 
(TTTft- fO]) reduces to 

ed t fi + 2pd t fi + ft d t p- Ap = infix(0,T), (1.11) 

6d t p-Ap + f'(p) = p J infix (0,T), (1.12) 

d vl i = d u p = on T x (0, T), (1.13) 

p(-,0) = p an d p( • , 0) = po in fi. (1-14) 

Note that h might attain its lower bound for some significant values of p, that is, for some 
p's lying in the domain of fa: actually, this was the case for h defined as in (II. 9p over 
the interval [0, 1], that is, over the effective domain of both potentials in (11.101) and (11.71) . 
We were prompted to generalize (11.91) as in (11.61) by an interesting remark of Alexander 
Mielke, when one of us was lecturing on our results, namely, that the behavior of 

h(p) = p + small parameter 

is different in a right neighbourhood of (h(p) m 0) than in a left neightbourhood of 1 
(h(p) « 1), whereas assuming only that h be bounded from below allows for many other 
instances like, e.g., a specular behavior of h around the extremal points of the domain of /. 

Returning now to (ll.ip - (ll.4l) . we set 

g(p) := h(p) - ~ > for all p G D(fa), 

and we reformulate our initial and boundary value problem as follows: 
to find p, p, and £, so as to solve 

(e + 2g(p))d t p + pg'(p)d t p-div(K(p)Vp) =0 infix (0,T), (1.15) 

6d tP -Ap + ^ + f^p)=pg'(p), with £G/3(p), infix(0,T), (1.16) 

(/t(p)Vp) • u\ T = d u p\ r = on T x (0, T), (1.17) 

p( • , 0) = po an d p( • , 0) = po in Q. (1-18) 

The global-in-time existence result we derive is more general than in [S], for two reasons: 
because it holds when the potential / includes a multivalued graph /3 (possibly with 
vertical segments, e.g., for fa as in (II. 7p ) and only exploits the monotonicity property of 



1 Note that, according to whether or not a\ > 2ct2, it turns out that / is convex in the whole of [0, 1] 
or it exhibits two wells with a local maximum at p = 1/2. 
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(3; and because the atom mobility n(p) is allowed to depend in a generic nonlinear way 
on the chemical potential. Note that problem fll.15p - fll.18p may become singular with 
respect to p due to the possible occurrence of singularities of /3; on the other hand, it may 
be also degenerate with respect to p since n(p) is allowed to vanish at p = 0. 

Our uniqueness result is also more general than in [8], because nonsmooth potentials 
could not be handled with the technique there used, consisting in testing the difference 
of two equations (11. lip by the time derivative of the difference of the two p components; 
however, just as in [8], the proof is achieved under the assumption that n(p) = a constant. 

Some directions for future research have already been explored by us under less general 
circumstances than those considered here: the long time behavior of system fll.lip - fll.14p 
and the structure of the relative omega-limit set have been analysed by us in [8] and in [9] , 
where we also dealt with the asymptotics of fll.lip - fll.14p as e — > and found a weaker 
solution in the singular limit. Moreover, in [10] and [12] we studied two optimal control 
problems for systems similar to (jl.llj) — (j!.14j) : a distributed control problem in JTU] and 
a boundary control problem in [12] . Finally, in [IT] we developed an existence theory for 
problem fll.ip - fll.4l) when atom mobility is allowed to depend on both p and p. 

This paper is organized as follows. In the next section, as anticipated, we discuss 
the physical features of the phase segregation model we adopt. In Section 3, we state 
our assumptions and results with the necessary mathematical accuracy; since here we do 
not take up asymptotic procedures, without loss of generality we set e — 1 in (11.151) and 
5 = 1 in (I1.16p . The existence of solutions to problem fll.15p - fll.18p is proved in Section HI 
their regularity properties in the successive section. Our last Section is devoted to the 
uniqueness proof. 



2 Short reasoned history of our mathematical model 

The nonstandard phase-field model fll.lip - fll.14p can be regarded as a variant of the classic 
Cahn-Hilliard system for diffusion-driven phase segregation by atom rearrangement: 

d tP -KAp = 0, p = -Ap + f(p). (2.1) 

Apart for the harmless choice n = 1 for the mobility modulus in (11.111) . one finds in 
fll.lip - fll.14p two awkward nonlinear terms involving time derivatives. Usually, equations 
(12.11) are combined in order to obtain the well-known Cahn-Hilliard equation: 

d t p = KA(-Ap + f'(p)). (2.2) 

Fried and Gurtin's generalization of Cahn-Hilliard equation. In [T5], [T7] , abroad 
generalization of (12. 2p was achieved by proposing the following: 



(i) to interpret the second of (12.11) as a balance of microforces: 

div£ + ii + 7 = 0, 



(2.3) 
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where the distance microforce per unit volume is split into an internal part it and 
an external part 7, and the contact microforce per unit area of a surface oriented 
by its normal n is measured by £ ■ n in terms of the microstress vector 

(ii) to regard the first of (12. II) as a balance law for the order parameter: 

d tP = - div h + a, (2.4) 
where the pair (h,a) is the inflow of p; 

(iii) to demand that the constitutive choices for tt, £, h, and the free energy density if;, 
be consistent in the sense of Coleman and Noll j5] with an ad hoc version of the 
Second Law of Continuum Thermodynamics: 

d t i[> + (tt - p)d tP - i ■ V{d tP ) + • Vp < 0, (2.5) 

that is, a postulated "dissipation inequality that accommodates diffusion" (cf. equa- 
tion (3.6) in [17]). 



In [17], the following set of constitutive prescriptions was shown to be consistent with (iii): 
{ i/j = $(p, Vp), 7?(p, Vp, fi)=fi- «9 p ^(p, Vp), |(p, Vp) = <9v>(p, Vp) } . (2.6) 
Moreover, it was presumed that 

h = -MVp, with M = M(p, Vp, p, Vp), (2.7) 
with the tensor-valued mobility mapping M satisfying the residual dissipation inequality 

Vp • M(p, Vp, p, Vp) Vp > 0. 

With the help of (12.31) . (12.41) . (I2.6p . and (I2.7p i . a general equation for diffusive phase 
segregation processes is arrived at, namely, 

dtp = div (MV (d${p, Vp) - div (<9 Vp ^(p, Vp)) - 7) ) + a. (2.8) 

The classic Cahn-Hilliard equation (I2.2p is obtained from (I2.8P by taking 

^(p,Vp) = /(p) + i|Vp| 2 , M = «i, (2.9) 

and by letting the external distance microforce 7 and the order-parameter source term a 
be identically null. 

An alternative generalization of Cahn-Hilliard equation. The Fried-Gurtin mod- 
el was well accepted in the mathematical community. In 2006, a largely modified version 
of it was proposed [19]: while the crucial step (i) was retained, both the order-parameter 

2 In [14], the microforce balance is stated under form of a principle of virtual powers for microscopic 
motions. 
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balance (12.41) and the dissipation inequality ( 12. 5 j) were dropped and replaced, respectively, 
by the microenergy balance 

d t e = e + w, e:=—dwh + a, w := —it d t p + £ ■ V(d t p), (2.10) 

and the microentropy imbalance 

d t r) > — div h + a, h := ph, a := p. (2-H) 

A further key feature of this new approach to modeling phase segregation by atomic 
rearrangement is that the microentropy inflow (h, a) is deemed proportional to the mi- 
croenergy inflow (h,a) through the chemical potential p, a positive field; consistently, the 
free energy is defined to be 

^:=e-n-\ (2.12) 

with chemical potential playing the same role as coldness in the deduction of the heat 
equation |f] 

Combining (l2~T0l -( l2~12|l yields 

dt^P < -Vdtil*- 1 ) +p- 1 h-Vfi-7rd tP + £- V(d t p), (2.13) 

an inequality that replaces (12. 5p in restricting a la Coleman-Noll the possible constitutive 
choices. On taking all of the constitutive mappings delivering tt, £, p, and h, dependent 
in principle on p, Vp, p, V/i, and on choosing 

V = Hp, Vp,p) = -p P + f(p) + i|Vp| 2 , (2.14) 
compatibility with ( 12.1 3D implies that we must have: 

{7r(p, Vp, fi) = -d p J/(p, Vp, p) = fi - f (p), 1 
?(p, Vp, p) = <9 Vp ^(p, Vp, p) = Vp, \ (2.15) 
V(p, Vp, p) = p 2 d^(p, Vp, p)= -p 2 p J 

together with 

/i(p, Vp, p, Vp) = H(p, Vp, p, Vp) Vp, Vp • i2"(p, Vp, p, Vp) Vp > 0. 

If we now choose for H the simplest expression H = nl , implying a constant and isotropic 
mobility, and if we once again assume that the external distance microforce 7 and the 
source a are null, then, with the use of ( 12.151) and ( 12.121) . the microforce balance ( 12. 3ft 
and the energy balance ( 12.101) become, respectively, 

Ap + p-/'(p) = (2.16) 

and 

2pd t p + pdtp - KAp = 0, (2.17) 

3 As much as absolute temperature is a macroscopic measure of microscopic agitation, its inverse - the 
coldness - measures microscopic quiet; likewise, as argued in [19], the chemical potential can be seen as 
a macroscopic measure of microscopic organization. 
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a nonlinear system for the unknowns p and p. 

Insertion of the parameters e and S. Let us compare systems (I2.16p -f l2.17p and 
(12. ip . Note that (I2.16P and f )2.ip 2 imply the same 'static' relation between p and p; 
instead. f)2.17p is rather different from ( 12.^ . for more than one reason: it is nonlinear; 
it features both time derivatives of p and p; and, in front of both d t p and d t p there 
are nonconstant factors that should remain nonnegative during the evolution. Thus, the 
system fl2.16p - fl2.17p deserves a careful analysis. 

We begun by attacking the problem as it was, prompted to optimism by the successful 
outcome of a previous joint research effort [6j [7] in which we tackled the system of Allen- 
Cahn type derived via the approach in [19] for no-diffusion phase-segregation processes. 
Unfortunately, the evolution problem ruled by ( I2.16p - fl2.17p turned out to be too diffi- 
cult for us. Therefore, we decided to study its regularized version f ll.lip - p.14p . where 
equations (II. lip and (I1.12p are obtained by introducing the extra terms e dtp and 5 dtp in 
(I2.17P and (I2.16p . respectively. Of course, the positive coefficients e and 5 were intended 
to be made smaller and smaller by way of an asymptotic procedure to be set up after the 
solvability of the regularized system were proved. 

Mathematically, the introduction of the e— term is motivated by the desire to have a 
strictly positive coefficient as a factor of d t p in (I2.17p . so as to guarantee the parabolic 
structure of equation (11. lip ; on the other hand, the introduction of the 5— term transforms 
(I2.16P into an Allen-Cahn equation with source p, and is strongly reminiscent of a sort 
of regularization already employed in various approaches to the so-called viscous Cahn- 
Hilliard equations (examples can be found in [2], [3], [161 HH1 ISHj and in the references 
therein) . 

It is also possible to make clear what additional physics the regularizing perturbations 
we introduced incorporate into the model. As to the term e d t p, it can be made to appear 
in the microenergy balance (11. lip by modifying as follows the choice for the free energy 



As to the term 5 dtp, it is enough to note that all is needed to make that term appear in 
the microforce balance (I1.12p is to add d t p to the list of state variables we considered to 
analyze the constitutive consequences of (12.131) . This measure brings in the dissipation 
mechanism typical of Allen-Cahn nondiffusional segregation processes, where dissipation 
depends essentially on (d t p) 2 , in addition to Cahn-Hilliard's |V/i| 2 — dissipation (cf. [T9~j); 
and it opens the way to splitting the distance microforce additively into an equilibrium 
and a nonequilibrium part, with 7r e9 = —dpijj(p,'Vp,p) = p — f'(p) the equilibrium part, 
just as in (12. 15[) x . and with n neq = —S d t p the nonequilibrium part. 



In this section, we state precisely the mathematical problem under investigation, fix our 
assumptions, and present our results. Let Q to be a bounded connected open set in M 3 
with smooth boundary T (the lower- dimensional cases can be treated with minor changes). 
We also introduce a final time T £ (0, +oo) and set Q := Q x (0,T). Moreover, we set 



in (J2HD: 




(2.18) 



3 Main results 
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for convenience: 

V:=H\Q), H:=L 2 {tt), and W := {v G H 2 (n) : d v v = on T}, (3.1) 

and we endow these spaces with their standard norms, for which we use a self-explanatory 
notation like || ■ \\y. For p G [1, +oo], we write || • || p for the usual norm in L P (Q). As no 
confusion can arise, the symbol || ■ || p is used for the norm in L P (Q) as well. Moreover, 
any of the above symbols for the norms is used even for any power of these spaces. We 
remark that the embeddings W C V C H are compact, since Q is bounded and smooth. 
As V is dense in H, we can identify H with a subspace of V* in the usual way (i.e., so as 
to have v *(u,v)v = (u,v)u for every u G H and v G V); the embedding H C V* is also 
compact. 

We are now concerned with the structural assumptions to set on our system. As the 
chemical potential is expected to be at least nonnegative, we assume that the function 
k is defined just for nonnegative arguments. However, one could study the more general 
mathematical problem of finding solutions whose component p might change its sign. In 
such has to be defined on the whole of R and must satisfy similar assumptions. 



We require that: 

k : [0, +oo) — > R is continuous, (3.2 

ft*, ft* G (0, +oo) and r* G [0, +oo), (3.3 

k(t) < ft* for every r > and ft(r) > ft* for every r > r*, (3.4 

K(r) := J Q r ft(s) ds for r > 0; K is strictly increasing, (3.5 

/ = /i + / 2 , /i :R^ [0,+oo], / 2 :R^R, p : R — )■ [0, +oo), (3.6; 

/i is convex, proper, l.s.c, and / 2 and are C 2 functions, (3.7 

/2, and g' are Lipschitz continuous, (3.8 

0:=dfi and vr := . (3.9 



In the following, D(j\) and D(/3) (C D(j\)) denote the effective domains of j\ and /?, 
respectively. 

Remark 3.1. We observe that our assumptions on / and ft allow for strong singularities 
(at the boundary of D(f3)) and a possible degeneracy (in a right neighbourhood of 0) in 
the equations for p and p, respectively. The former fact is clear. As far as the latter is 
concerned, we note that (13.51) is satisfied if and only if k is nonnegative and the set where 
ft vanishes has empty interior. So, r* = means uniform parabolicity for equation (11.11) . 
On the contrary, if r* > 0, the equation can degenerate, e.g., at the origin (or even in 
rather big set of small values). An example is given by ft(r) = tanhr m ~ 1 with m > 1. In 
this case, (II. ip roughly behaves like the porous medium equation (slow diffusion) in the 
region where p is small. 

Remark 3.2. It is known that any proper, convex and lower semicontinuous function 
is bounded from below by an afline function (see, e.g., [H Prop. 2.1, p. 51]). Hence, 
our assumption fx > looks reasonable, because one can suitably modify the smooth 
perturbation f 2 by adding a straight line. Moreover, the other positivity condition, g > 0, 
is just needed on the set D((3), while g can take negative values outside of D(j3). Finally, 
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(13. 8p implies that, within the range of relevant values of r, the functions f^r), g{r), and 
g'(r) grow at most linearly with respect to \r\, while /2( r ) grows at most quadratically 
in \r\. 

For the initial data, we postulate: 

p eV, po e W, p > and p G D((3) a.e. in Q, (3.10) 
and there exists some £o G H such that £o £ /3(po) a.e. in f2. (3-H) 

Since f\ is convex and is smooth, the above assumptions entail /(po) £ L l (^L). 

Now, we introduce the a priori regularity that we require from any solution (p,p, £) 
to our problem. Note that equation (jl.2p reduces for any given p to a rather standard 
phase- field equation. Therefore, it is natural to look for pairs (p, £) that satisfy 



p G W^O, T; #) n H\0, T; V) n L°°(0, T; W), (3.12) 

£gL°°(0,T;#), (3.13) 

and to solve the subproblem in a strong form, namely 

d tP - Ap + i + 7r(p) = p </(p) and £ G /3(p) a.e. in Q, (3.14) 

p(0) = po a.e. in f2. (3.15) 



We note that (13.121) also contains the Neumann boundary condition for p (see (13. ip for 
the definition of W). On the contrary, the situation is different for the component p. In 
the case of uniform parabolicity, i.e., if r* = 0, the coefficient ft(p) is bounded away from 
zero, and we can require that 

p G -f^ 1 (0, T; H) n L°°(0, T; V), p>0 a.e. in Q, (3.16) 
div(/<p)Vp) G L 2 (0, T; H), (3.17) 

and that p satisfy 

/ (l + 2g(p(t)))d t /M(t)v+ [ fi{t)g'{p{t))d t p(t)v+ [ «(p(t))Vp(t) • Vv = 
</n isi 

for every v G V and for a. a. t G (0,T), (3.18) 
p(0) = po a.e. in Q. (3.19) 

Thus, the equation holds in a strong sense, i.e., 

(1 + 2g(p)) d t n + fj.g'(p) d tP - div (/t(p)Vp) = a.e. in Q, (3.20) 

while the Neumann boundary condition is understood in the usual weak sense. Further- 
more, we observe that ( I3.16l) -( l3.18l) imply further regularity for p whenever k is smoother, 
thanks to the regularity theory of quasilinear elliptic equations. 

If instead we allow r* to be positive, such a formulation is too strong, since no sufficient 
information on the gradient Vp can be obtained. As a consequence, the same happens 
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for the time derivative d t p. For that reason, we rewrite equation (13.201) in the different 
form 

dt (1 + 2g{p)p) - pg\p)d t p - AK(p) = 0. (3.21) 

More precisely, we also account for the initial and Neumann boundary conditions and 
accordingly rewrite fl3.18p - fl3.19p . In conclusion, we require the lower regularity 

p G L°°(0, T; H), p > a.e. in Q, K(p) G H 1 ^, T; H) D L°°(0, T; V), (3.22) 

(1 + 2^6^(0^;^), (3-23) 
and replace fl3TT8|) - fl3TT9|) by 

(d t ((l + 2g(p))p)(t),v)- [ (pg'(p)d t p)(t)v + [ VK(p(t)) ■ Vv = 

Jn Jn 

for every v G V and for a. a. t G (0,T), (3.24) 

((1 + 2g{p))p) (0) = (1 + 2g(po)) p Q . (3.25) 

In this situation, f)3.2ip is satisfied in the sense of distributions, only. 

Remark 3.3. We observe that even the middle term of (13.241) is meaningful, as we 
immediately see. First, we note that 

peC°([0,T];C (n)) = C (Q), (3.26) 

directly from (13.121) and the compact embedding W C C°(f2) (see, e.g., [211 Sect. 8, 
Cor. 4]), whence g'(p) G C°(Q). Next, fl3T22|) and the embedding V C L 4 (fi) imply that 
K(p) G L°°(0,T;L 4 (fi)), whence also p G L°°(0, T; L 4 (fi)), since £T(r) behaves like r for 
big \r\ by (J32D- Finally, d t p G L°°(0, T; F). Therefore, pg'(p)d t p G L°°(0, T; L 4 / 3 (fi)). 
On the other hand, ?; G L 4 (Q) whenever v G V. 

Remark 3.4. Note that f l3T25|) makes sense because (1 + 2g{p))p G C°([0,T]; V*) (due 
to (I3.23P ). However, by accounting for f )3.12p and the regularity of g, we see that (13. 25 j) 
can be read in the simpler form (I3.19P also in this case. Indeed, the function (1 + 2g(p))p 
also belongs to L°°(0,T; H). As is well known (and easy to prove), this implies that it 
actually is an if -valued function which is weakly continuous, in addition. It easily follows 
that p enjoys the same property. 

Here are our results. The first two state that the strong and weak formulations are 
equivalent in the case = and that there exists a weak solution in the general case. 
Due to the former, the latter also proves the existence of a strong solution if = 0. Both 
results will be proved in Section |H 

Proposition 3.5. Assume f l3T2|) -f l3T9|) . f l3TT0|) - fl3TL|) . and r* = 0. Then, any triplet 
(p,p,0 satisfing <KY2\f - fl3~T5j) . f l3T2^|) -f lX25|) and solving problem (1334 ]) -(13TT5 ]) , fl£2§- 
Q3T25]) also satisfies (j3TTo]l - (J3TT5]) . 



Theorem 3.6. Assume (13. 2p - (13.91) and (13. 10p - (13.111) . Then, there exists at least one 
triplet (n,p,£) satisfing fl3TT2D - fl3TT5|) . fl322D -(JS23D and solving problem fl3TT4D - fl3TT5]) . 

(JXMD-(JS25D. 
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We notice that no further assumptions are needed to ensure boundedness for p, due 
to f)3.26p . As far as the first component is concerned, we have the following boundedness 
result. 

Theorem 3.7. Assume (13T2| -(1331) . (13~T0l) -l ETTTj) . and let p G L°°(n). Then, the com- 
ponent p of any triplet (/i,p, £) satisfing ( 13. 12ft -( 13. 13ft . ( 13.22ft -( 13. 23ft and solving prob- 
lem (I5TTD - (I3TT5]) and (I3^4j) - (I3T25D is essentially bounded. 

The next results hold if we assume that k is constant. We notice that we could weaken 
this assumption in our regularity result, while we are not able to do the same as far as 
uniqueness is concerned, unfortunately. In order to simplify the regularity proof, we take 
K = 1 at once. In the forthcoming Remark 15.41 we will sketch how to deduce even further 
regularity. 

Theorem 3.8. Assume (13T2D -(I3T9D . (jSTTUft - ({5TTT|) , p G W, and k = 1. Then, any triplet 
(/i,p,0 satisfing (I3TT2|) - (13331) . (I336|) and solving problem (13?L4j) - (13351) and (I3TT8D - (I339|) 

enjoys the regularity property 

p G W 1,p (0, T; H) fl L p (0, T; W) for every pG [1, +oo). (3.27) 

Theorem 3.9. Assume (I3T2|) - (I3T9|) . (l3TT0]l - (j3TTTl ). p G W, andn=l. Then, the triplet 
(fi,p,0 satisfing (|3TT2]1 - fl3TTB1) . <KW\i and solving problem (1334]) - (13TT51) and (133B1 - (13391 
is unique. 

Throughout the paper, we account for the well-known embedding V C L P (Q) for 
1 < p < 6 and the related Sobolev inequality: 

II v ||p < Clkllv for every v £ V and 1 < p < 6, (3.28) 

where C depends on Q only. Moreover, we recall that the embeddings V C L 4 (f2) (more 
generally V C L P (Q) with p < 6) and W C C (fi) are compact and use the corresponding 
inequality 

IMU < e||Vu + C e ||?;||# for every v G V and £ > 0, (3.29) 
where C e depends on Q and e, only. Furthermore, we make repeated use of the notation 

Q t := tt x (0, t) for t G [0, T], (3.30) 

and of the well-known Holder inequality and the elementary Young inequality 

ab < ea 2 + — b 2 for every a, b > and e > 0. (3.31) 

Finally, throughout the paper we use a small-case italic c for different constants that may 
only depend on Q, the final time T, the shape of the nonlinearities / and g, and the 
properties of the data involved in the statements at hand; a notation like c £ signals a 
constant that depends also on the parameter e. The reader should keep in mind that 
the meaning of c and c £ might change from line to line and even in the same chain of 
inequalities, whereas those constants that we need to refer to are always denoted by capital 
letters, just like C in fl3T2g|) . 
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4 Existence 

In this section, we first show the equivalence result stated in Proposition 13. 51 Then, we 
prove Theorem 13. 6| which ensures the existence of a weak solution. 

Proof of Proposition Iff. 51 As r* = 0, we have k(t) > «;* for every r > 0. This implies 
that the inverse function K~ x : [0, +oo) — > [0, +oo) is Lipschitz continuous. Hence, (13 .22 j) 
implies that 

p = K~\K{n)) G H l (0, T- H) n L°°(0,T; V), 
i.e., that (13.161) holds. In particular, we can write 

VK(ji) = K(ji)Vfi and d t ((l + 2g{p))p) = pd t (l + 2g(p)) + (1 + 2g(p))d t p 

and thus replace the weak formulation by the strong one. Next, we note that (I3.18P implies 
that (13.201) holds in the sense of distribution, whence (13.171) follows by comparison. Finally, 
(I3.19P holds even in the general case, as we have observed in Remark 13.41 □ 

Now we prove Theorem 13.61 Even though our proof closely follows the argumentation 
of [8] , we present the whole procedure and sometimes give some detail, since the changes 
with respect to the quoted paper are spread in the whole calculation. The starting point 
is an approximating problem which is still based on a time delay in the right-hand side 
of (13341) . Namely, we define the translation operator T r : L\Q,T;H) -> L\Q,T;H) 
depending on a time step r > by setting, for v G L 1 (0,T; H) and for a.a. t G (0, T), 

(7 T v)(t) := v(t - t) if t > r and (7 T v)(t) := p if t < r (4.1) 

(but the same notation 7 T v will be used even for functions v that are defined in some 
subinterval [0,T'] of [0,T]), and replace \i by T r /i in (I3.14p . essentially. However, we 
modify the equation for \i at the same time. Precisely, we force uniform parabolicity and 
allow the solution to take negative values, if possible. To do that, we define k t : R — > R 
and the related function K T used later on by 

K T (r) := Ac(|r|) + r and K T (r) := / k t (s) ds for r G R. (4.2) 

Jo 

So, the approximating problem consists of the following equations 

(l + 2g(p T )) dtfir + ^rg'(Pr) d t p r - div (k t (fi T ) V/i r ) = a.e. in Q, (4.3) 
d t pr - Ap r + £ r + n(p T ) = (7 T p T ) g'(p T ) and £ r G /3(p T ) a.e. in Q, (4.4) 

complemented with the homogeneous Neumann boundary conditions for both p T and p T 
and the initial conditions p T (0) = po and p T (0) = po- For convenience, we allow r to take 
just discrete values, namely, r = T/N, where is any positive integer. □ 

Lemma 4.1. The approximating problem has a solution (/i T , p r , £ r ) satisfying the ana- 
logues of fl3TT2|) - fl3TTk and fl3TT6l) - fl3T?]) . 

Proof. We just give a sketch. As in [S], we inductively solve iV problems on the time 
intervals I n = [0,t n ], n = 1, . . ., N, by constructing the solution directly on the whole 
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of I n at each step. Namely, given p n -i, which is defined in f2 x J n _i, we note that 7 T p n _i 
is well defined and known in Q x I n (even in the starting case n — 1) and solves the 
boundary value problem for p n given by the phase-field equations 

d t p n - Ap n + £ n + vr(p n ) = (J T p n ) g'(p n ) and G /3(p n ) in x J n , (4.5) 

complemented with the boundary and initial conditions just mentioned for p T . Such a 
problem is quite standard and has a unique solution in a proper (rather weak) functional 
analytic framework. Once p n is constructed, we solve the parabolic equation 

(l + 2g(p n )) d t pn + fi n g'(p n ) d t p n - div (k t (p n ) V p n ) =0 in Q x I n , (4.6) 

together with the boundary and initial conditions prescribed for p T . We note that g > 
and K T (r) > r for every r G M, so that the equation is uniformly parabolic. Therefore, 
the problem to be solved has a unique solution in a proper space provided that the 
coefficient g'(p n )d t p n is not too irregular. So, we should prove that, step by step, we get 
the right regularity for p n and p n . This could be done by induction, as in [5], with some 
modifications due to our more general framework. We omit this detail and just observe 
that the needed a priori estimates are close (and even simpler, since r is fixed here) to the 
ones performed later on in order to let r tend to zero. The final point is p n > 0. We give 
the proof in detail. We multiply equation (14. 6[) by — p~ := — (— p n ) + , the negative part 
of p n , and integrate over Q t with any t G I n . We observe that 

[(l + 2g(p n (t)))d t p n + p n g(p n )d t p n ] {-p~) = )^d t ({l + 2g{p n )) \p~\ 2 ) . 
Hence, by using p > 0, and owing to the boundary condition, we have 

~ / (1 + 2g(p n (t))) \p-{t)\ 2 + / K T {p n )\Vp-\ 2 = for every t G I n . 
A Jn JQ t 

Since g and k t are nonnegative, this implies p~ = 0, that is, p n > a.e. in Q x I n . Once 
all this is checked, the finite sequence (/i n , p n , £„), n = 1, . . . , N, is actually constructed, 
and it is clear that a solution to the approximating problem we are looking for is obained 
by simply taking n = N. □ 

Although the solution to the approximating problem is unique, we do not need unique- 
ness in the following and just fix a solution (p T ,p T ,^ T ) for each r. Our aim is to let r 
tend to zero in order to obtain a solution as stated in Theorem 13.61 Our proof uses com- 
pactness arguments and thus relies on a number of uniform (with respect to r) a priori 
estimates. Clearly, in performing them, we can take r as small as we desire, and it will 
be suitable to assume that r < k*. In order to make the formulas more readable, we 
shall omit the index r in the calculations, waiting for writing p T and p T only when each 
estimate is established. 

First a priori estimate. Let us test H4.3[) by p T and point out that 



[(1 + 2g{p)) d t p + p g'{p) d t p] p = -d t [(l + 2g( P ))p 2 ] . 
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Therefore, by integrating over (0, t), where t e [0,T] is arbitrary, we obtain 
(l + 2^(p(t)))|/i(t)| 2 + / Kr(ji)\Vii\ 2 = f(l + 2g{po))f4. 



Hence, we recall that g > and observe that K 2 (r) < 2k*k t (t) for every r e R by (13 .4p 
and r < k*. We conclude that 

\\Pt\\l°°(0,T;H) + ||^t(/^t)|U 2 (0,T;V) < C. (4.7) 

Actually, we have proved more, namely 

|| JC(/ir) || l2(o,t ; v) < c where K* T (r) := f (k t (s)) 1/2 ds for r G R. 

Jo 

Moreover, we observe that K has a linear growth, so that (14.71) also yields 

\\K T (p T )\\ L °°(0,T;H) < c. (4.8) 

An implication of (H77|l -(Oj). along with (j4TTJ and (ETT0I) . is 

||T T ^ T ||£oo(0,T;H) + ||'J'r-ft'r(A i r)IU° o (0,T; J H')nL 2 (0,T;y) < C. (4.9) 

Consequence. The Sobolev inequality (I3.28P and estimate (14. 7p imply that 

||JC(A t r)||L 2 (0,r ; L6(f7)) < C. 

On the other hand, (13.41) implies that K T (r) > — c for every r > 0. We deduce that 

||A*r|U2(0,T;i6(n)) < c (4.10) 

Second a priori estimate. Let us add p T on both sides of (14.41) and test by d t p T . We 
have that 

\d tP \ 2 + l\\p(t)\\v+ f fMt)) 

= ^IIPo|| 2 y+ / /(Pb) + i / (p 2 (t)-2/ 2 (p(t))) + / g'(p)(7 T p)d tP 
L Jn L Jn JQ t 

<c + c [ \p{t)\ 2 + - / |^p| 2 + cp T p,f Loo ^ T . H) , 
Jn 4 JQt 

for every t G [0,T]. In view of the chain rule and Young's inequality (13.31 j) . we observe 
that 

c j \ P (t)\ 2 <c [ |p | 2 + i I \d t p\ 2 + c f\\ P (s)\\lds. 
Jn Jn 4 JQ t Jo 

Hence, as fi is nonnegative, on account of (14. 9p . and with the help of the Gronwall lemma, 
we deduce that 

\dtp\ 2 + \\p(t)\\v+ I fMt))<c 
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Therefore, we obtain: 

\\Pr\\m(0,T;H)nL°°(Q,T;V) < C and || f(p T ) || L oo (0jT;L i {n)) < c. (4.11) 

Third a priori estimate. We rewrite (14.41) as 

- Ap + p(p) 3 -d t p - 7r(p) + {%p)g'{p) 

and note that the right-hand side is bounded in L 2 (0,T; H), thanks to the Lipschitz 

continuity of ir and g' and to the previous estimates. By a standard argument (formally 
test by either — Ap or /3(p) and use the regularity theory for elliptic equations), we first 
recover that 

IIM«)IIIr + IKOOII! < 2||-M«) - *(pOO) + ((7^)g\p))(s)\\l (4.12) 
for a. a. s G (0,T), and finally conclude that 

\\Pt\\lz(o,T;W) < c and ||£ T \\l*(o,t;H) < c. (4.13) 

Fourth a priori estimate. Our aim is to improve the estimates (14. lip and ( 14.131) . 
To this end, we proceed formally, at least at the beginning, for the sake of simplicity. 
However, this procedure could be made rigorous by suitably regularizing equation (I4.4p 
(with respect to p, only, i.e., keeping p fixed), the main tool being the Yosida regularization 
of maximal monotone operators (see, e.g., [U p. 28]; see also the proof of Lemma 3.1 of [5] 
for a further regularization). Such a theory yields, in particular, the estimate 

||3tu(0)||fr< U(0) + Ap || H + mm\\r]\\ H (4.14) 

for the unique solution (u, u) to the equations 

d t u — Au + u = ip := g'(p)7 T p — 7r(p) and uj G (3(u), 

complemented with the same initial and boundary conditions as those prescribed for p. 
Note that in (14. 14ft /3 is understood as the induced maximal monotone operator from H 
to H. Observe that (u,oj) = (/?,£); then the application of (I4.14p . in combination with 
our assumptions on p (see (13.111) . in particular), leads to 

\\dtp T (0)\\ H < c(\\ho\\h + \\po\\w + 1 + Uo\\h)= c. (4.15) 

We use (14.151) in the subsequent calculation, where we proceed formally, as announced 
(however, our procedure becomes completely rigorous after a while). In particular, we 
write 0(p) in place of £ and treat (3 as if it were a smooth function. We differentiate ( I4.4p 
with respect to time and test the resulting equation by d t p. We obtain: 



\ [ \dtP(t)\ 2 + ( |V«9 t p| 2 + / (3'(p)\d tP \ 2 

A JQ JQt JQt 

\(d tP )(0)\ 2 - [ ApMp\ 2 + [ g"{ P ){%p)\d t p\ 2 



+ / g\p)d t {7 T p)d t p 
JQt 

<\f \(dtp)(0)\ 2 + c / (1+T r p)|^p| 2 + / g'(p)d t (7 T p)d t p. (4.16) 
z Jn JQ t JQt 
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We treat each term on the right-hand side separately. The first one is estimated by ( 14.151) . 
In order to deal with the second one, we first account for the Holder inequality. Then, 
we also invoke (14. 7p . the compact embedding V C L 4 (f2) (see (I3.29P ). and (14. lip . For 
every e > we infer that 

/ (l+7rri\d tP \ 2 < r||l + (^)(s)|k||^)lll^ 
JQt JO 

< c [ \\d tP (s)\\lds < e [ \Vd t p\ 2 + c £ [ \d tP \ 2 

JO JQt JQt 

<e [ \Vd t p\ 2 + c £ . (4.17) 

JQt 

The estimate of the last term of (14.161) needs much more work. We recall that 7 T p is 
constant with respect to time on the interval (0, r) and first compute d t p from (14. 3p . Then 
we integrate by parts and repeatedly exploit the Holder, Sobolev, and Young inequalities. 
We obtain: 

g\p)d t (7 T p) d t p = / / d t fi(s)g'(p(s + T))d t p(s + T)ds 



Q t Jo Jn 

t-T 



o J n l + 2g(p(s)) 



[div(K T (p)(s)Vp(s)) - fi(s)g'(p(s))dtp(s)]d t g(p(s + T))ds 



t-T 



, w , s „d t q(p(s + T)) , 
l + 2g(p(s)) 



o Jn 



-n '™s*.» T)) *™'™- +r) * (4 - is) 

and now treat the last two integrals separately, by accounting for our structural assump- 
tions. We have 

Jo Jn 1 + 2#(p(s)) 

Jo Jn l + 2g{p(s)) 

<c f [ \VK T (p)(s)\ \Vd t p(s + r)\ds 
Jo Jn 

+ cf [ \VK T (p)(s)\ \Vp(s)\ \d t p(s + r)\ds 
Jo Jn 

+ c [ 7 \VK T (p)(s)\ \Vp(s + r)| \d t p(s + t)\ ds. (4.19) 
Jo Jn 

For every e G (0, 1) there holds 

-t-T 



VK T (p)(s)\ \Vd t p(s + T)\ds 

\Vd tP \ 2 + c £ [ \VK T (p)\ 2 <el \Vd t p\ 2 + c £J (4.20) 
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thanks to (j4.7|) . On the other hand, we also have 

f 7 |Vif T 0*)(a)||Vp( a )||^( a + r)|d« 

< / T \\VKM(s)\\ 2 \\Vp(s)\\4d t p(s + r)\\ 4 ds 
Jo 

<e f\\d t p{s)\\ 2 v ds + c £ f T \\VKM(s)\\ 2 H \\V P (s)\\ 2 v ds 
Jo Jo 



Qt 



<e \Vd t p\ 2 + c \d t p\ 2 + c £ \\VK T (p)(s)\\ 2 H \\Vp(s)\\ 2 v ds 



t-T 



<e f \Vd tP \ 2 + c + c £ [ r ||V^(//)(s)||^||Vp( S )||^ S . (4.21) 

jQt Jo 

In the last inequality we have used (14.111) . However, the above estimate has to be im- 
proved. To this end, we use the regularity theory for linear elliptic equations and estimates 
(14.91) and (14. lip again. Indeed, with the help of (I4.12p we have: 

\\Vp(s)\\ 2 v < c{\\p(s)\\ 2 v + ||M*)HSr) < c{\\d t p(s)\\ 2 H + 1). 
Therefore, the above estimate becomes 

-t-T f 

\WK T (p)(s)\\Vp(s)\\d t p(s + r)\ds 




'o Jn 



<e f \Vd t p\ 2 + c £ f \\VK T (p)(s)\\%\\d t p(s)\\%ds + c E . (4.22) 

JQt Jo 

Analogously, one shows that 

/ 7 \VK T (p)(s)\\Vp(s + r)\\d t p(s + r)\ds 
Jo Jn 

<e [ \Vd tP \ 2 + c £ ! \\V(7 T K T (p))(s)\\ 2 H \\d t p( S )\\ 2 H ds + c E . (4.23) 
JQt Jo 

Thus, by collecting (|420|) and rt4T22]) - d423|) . we deduce that (T4TT9D yields 

f 7 K T (p)(s)Vp(s) • V 9tP< f + T) ds < e j \Vd tP \ 2 
Jo Jn l + 2g[p(s)) J Qt 

+ c £ jf (\\VK T (p)(s)\\ 2 H + \M7 T K T (p))(s)\\ 2 H ) \\d tP (s)\\ 2 H ds + c £ (4.24) 

for every e > 0. Let us come to the last term of ( 14.181) . By using the compacness 
inequality (13.291) . and (14.71) as well, we have 

r g\p{s))g'{p{s + r)) 

p{s)d t p{s)d t p{s + r) ds 




o J n l + 2g(p(s)) 



t-T 



<C \\p(s)Ud t p(s + T)\\ 4 \\d t p(s)\\ 2 ds 

Jo 

<e f T \\ dt p(s + r)\\ 2 v ds + c £ f\\ti*)\\l\\0trt*)\\%d8 
Jo Jo 

<e I \\Vd t p{s)\\ 2 H ds + c + c £ I \\p(s)\\l\\d t p(s)\\ 2 H ds. (4.25) 
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Therefore, due to g^D and fTOBD . fTCTS]) becomes 

g'(p)d t (7 T fi)d t p<2e [ \Vd tP \ 2 + c £ [ \\pi(s)\\l\\d t p(s)\\ 2 H ds 

c £ [ f||Vir r (/i)( S )||i+||V(T r AV(/i))( S )||^||^( S )||^rf S + c £ 



(4.26) 



At this point, we combine (14.151) . (I4.17p . (14.261) with (14.161) and choose e small enough. 
As the last integral on the left-hand side is nonnegative since f\ is convex, we obtain 

\dtp(t)\ 2 + [ \Vd tP \ 2 <c f <P(s)\\d tP (s)\\ 2 H ds + c 

Q JQt JO 

where <f>(s) := \\f,(s)\\ 2 + \\VK T (p)(s)\\ 2 H + \\V(7 T KM)(s)\\ 2 H . 
As G -^ 1 (0, T) by (14. 7p and (I4.10p . we can apply the Gronwall lemma and conclude that 

||^PT||L°°(0,T;H)nL 2 (0,T;\/) < c - (4.27) 

Consequence. We have — Ap T + £ r = ip := — d t p T + g'(p T )7 T p T G L°°(0,T; H) due to 
(14.71) and (14.271) . Therefore, by a standard argument (formally multiply by —Ap T at any 
fixed time), we deduce that both —Ap T and £ r belong to L°°(0, T; H). Owing to elliptic 
regularity, we conclude that 



\Pt\\l°°{o,t-,w) < c and \\£t\\l°°(o,T;H) < c, (4.28) 



whence also 



||Pt|U-(Q) + \\9(Pt))\\l°°(Q) + W(Pt))\\l°°(Q) + h{pr))\\ L ^(Q) < c, (4.29) 

due to the continuous embedding W C L°°(Q) and the continuity of g, g', and ir (note 
however that g' is a bounded function since g is Lipschitz continuous). 

Fifth a priori estimate. We write (14. 3 p as d t ((l + 2g(p))fx) = AK T (p) + g'(p)fid t p. 
Thus, we have for every v G L 2 (0, T; V) 



Q 



d t ((l + 2g(p))p)v 



Q 



VK T {fj,)-Vv+ / g\p)pd t pv 
Jq 



< (\\Kt(h)\\l2(0,T;V) + c\\dtp\\L°°(0,T;H) \ \p\\L 2 (0,T;L 4 (Q))) \\ V \\ L 2 {0,T;V) ■ 



By accounting for (14TT1) . (H~27]) . and (T4TTUT) . we deduce that 

\\d t ((l + 2g(p T ))fi T ) \\mo,T;V*) < c 

Sixth a priori estimate. We test (14.31) by d t K T (p) = n T (p)d t p and obtain 

1 



(4.30) 



{l + 2g{p))K T {p)\d t p,\ 



|VA- T (/x )P 



\VK T {ti{t))\" 



g'(p)dtpnd t K T (ti) 



(4.31) 
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for every t G (0, T). Note that the first term on the left-hand side can be estimated from 
below as follows 

/ (l + 2g(p))KM\d t p\ 2 > [ ^4\d t p\ 2 = ^[ \d t KM\ 2 - (4.32) 
JQt JQt ZK Zk JQt 

Now, we deal with the right-hand side of (I4.3ip . The first term being trivial thanks to 
()3.10P i . we come to the second one. We invoke the Young, Holder, and Sobolev inequalities 
and have 

- / g'(p)d t ppd t K T (p)<-L [ \d t KM\ 2 + c [ \\fi(s)\\l\\d t p(s)\\lds 
JQt Ak JQt JO 

<-L/ \d t KM\ 2 + c [)\p(s)\\l\\d tP (s)\\ 2 v ds. (4.33) 

4K JQt Jo 

Now, we observe that (13 .4p yields K T (r) > K*r — c* for every r > 0, where c* depends on 
the structural assumptions, only. Hence, by owing to (14.81) as well, we deduce 

\\p(s)\\l < c(\\K T (p)(s)\\l + 1) < c(\\K T (p)(s)\\ 2 v + 1) 

< c\\VK T (p)(s)\\ 2 H + c\\KM(s)\\ 2 H + c < c\\VK T {^)f H + c 

for a.a. s G (0,T). By combining (l4T3"2l and fl4T3"31 with (OB . we obtain 

f \d t KM\ 2 + 1 - I \VKMt))\ 2 <c + c f<f,(a){\\VK T (jM)(8)\\ 3 H + l)da 

where 0(s) := ||9 t p(s)||y. 
As G L 1 (0,T) by (14.271) . we can apply the Gronwall lemma and conclude that 

\\Kr(pT)\\m(0,T;H)nL^(0,T;V) < C. (4.34) 

Consequence. By arguing as we did for (I4.10p . we derive that 

||Mr|U°°(0,T;L6(n)) < C. (4.35) 

Limit and conclusion. By the above estimates, there exist a triplet (p,p, £), with 
p > a.e. in Q, and functions k and £ such that 



Pt P 


weakly star 


in 


L~(0,T;L 6 (fi)), 


(4.36) 


Pr^P 


weakly star 


in 


L°°(0,T; H/), 


(4.37) 


d t pr -> <%p 


weakly star 


in 


l°°(o,T; F)nL 2 (o,r ; y), 


(4.38) 




weakly star 


in 


L°°(0,T;tf), 


(4.39) 


K T (p, T ) —7- fc 


weakly star 


in 


^(0,7; H)nL°°(0,T; V), 


(4.40) 


C T := (1 + 2g(p T )) fJ , T -> C 


weakly star 


in 


H^O, T; V*) n L°°(0, T; L 6 (fi)), 


(4.41) 



at least for a susequence r = Tj\0. By (l4.37p - (l4.38p . (14 .40 p . and the compact embeddings 
W C C°(Q) and V C if, we can apply well-known strong compactness results (see, e.g., 
[211 Sect. 8, Cor. 4]) and have that 

p T -)■ p strongly in C°(Q) (4.42) 

K T (p T ) k strongly in C°([0, T}; H) and a.e. in Q. (4.43) 
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The weak convergence (14.391) and (14.421) imply that £ G /3(p) a.e. in Q, as is well known 
(see, e.g., [H Prop. 2.5, p. 27]). The strong convergence (I4.42p also implies the Cauchy 
condition (I3.15P and that <fr(p T ) —> <f)(p) strongly in C°(Q) for every continuous function 
<fi : K — y R. We can apply this fact to the functions g, g', and tt (see (I3.8P ). In particular, 
we infer that p T g'(p T ) has some weak limit in L°°(0, T; L 6 (f2)): thus, we can identify it 
with the help of (I4.36P and conclude that (I3.14p holds. Now, we prove that p T converges 
to p a.e. in Q. To this aim, we note that K~ x converges to K~ l uniformly on [0, R] for 
every R > 0. Hence, we see that (I4.43P implies p T — > K~ l {k) a.e. in Q. By comparison 
with (14.361) . we deduce that p T — > p a.e. in Q. Next, let us deal with the subproblem 
for p. The identity K^ 1 (k) = p just proved means that k = K(p). From the convergence 
almost everywhere of p T to p we also infer that ( T converges to (1 + 2g(p))p a.e. in Q, 
whence ( = (1 + 2g(p))p by comparing with (14.41 j) . The last term to be identified is the 
limit of r] T := p T g' \p T )dtp T - Precisely, we prove that r\ T converges to r\ := pg'(p)d t p weakly 
in some L p -type space. We observe that (I4.36P and the convergence almost everywhere of 



p T imply that 

p T —7- p strongly in L p (0, T; L q (Q)) for every p < +oo and q < 6 (4.44) 

as is well-known (via the Severini-Egorov theorem). By choosing, e.g., p = q = 4 and 
combining with the weak star convergence of dtp T in L°°(0, T; H) (see (I4.38P ) and the uni- 
form convergence of g'(p T ), we deduce that r] T converges to r] weakly in L 4 (0,T; L 4 / 3 (f2)), 
thus weakly in L 2 (0, T; L 4 / 3 (f2)). At this point, it is straightforward to derive (13.241) in 
an integrated form, namely 

(d t ((l + 2g(p))p)(t),v(t))dt- [ pg'(p)d t pv+ [ VK(p)-W = (4.45) 

J Q JQ 

for any v G L 2 (0,T; V) C L 2 (0,T; L 4 (f2)), whence also the time-pointwise version (I3.24[) 
itself. Finally, (1PTD implies that ( T ( strongly in C°([0, T]; V*), thus, Cr(0) -> C(0) 
strongly in V*, so that the Cauchy condition (13.251) is verified as well. This concludes the 
proof. 

5 Further properties 

In this section, we prove Theorems 13.71 and 13.81 and make some remarks on the regularity 
of solutions. As far as the first result is concerned, we adapt the argument used in [8]. 
However, as the first estimate of the proof has to be derived in a different way, we prepare 
a technical lemma. 

Lemma 5.1. Assume 

u G i/ 1 (0, T; V*) fl L°°(0, T; H) and u + G -£^(0, T; H) R L°°(0, T; V), 
and let 7 G H l (0, T; V) n L°°(0, T; W). Then, for every t G [0, T] , we have that 

(d t u(s),~/(s)u + (s)) ds = / wyd t u + . (5.1) 
JQ t 
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Proof. As in Remark 13. 44 u is a weakly continuous if- valued function and the pointwise 
values of u and u + make sense. We start from the formula 

(d t u(s), v(s)) ds = (u(t), v(t)) - (u(0), v(0)) - / (u(s), d t v(s)) ds 

Jo 

u(t)v(t)- I u(0)v(0)- / ud t v, 
Jn JQ t 

which holds if v G if 1 (0,T; V). By an easy regularization, one sees that it still holds if 
v G i/ 1 (0, T; H) fl L 2 (0,T; V). Now, by applying our assumptions on u + and 7 and also 
owing to the Sobolev inequality (I3.28p . we have 

u + G L°°(0, T; L\Q)), V« + G L 2 (Q), d t u + G L 2 {Q) 

7GL°°(g), V 7 6r(0,T;L 4 (n)), %GL 2 (0,T;I 4 (fi)). (5.2) 

It follows that 7U+ G H^O^T; H) n £ 2 (0,T; V). Therefore, we obtain 

{d t u(s),'j(s)u + (s)) ds 



u{t)^{t)u + {t) - / m(0)7(0)m + (0) - / m^(7m + ) 
M(t)7(t)w + (t) - / m(0)7(0)m + (0) - / n + (n + % + 7^n H 



7 (t)|u+(t)| 2 - / 7 (0)|u+(0)| 2 - / \ u +\ 2 d tl - 1 d t \u+\ 2 + u + 7 «9 t u + 
n Jn JQ t JQt JQt 

wyd t u + 

Qt 

for all t G [0,T], and the lemma is proved. □ 

Proof of Theorem \3. 1\ Set /!q := max {1, ||wo||oo}- I n the proof performed in [8] the 
quantity J n | (//(£) — k) + \ 2 + f Qt | V(/i — k) + \ 2 is estimated for any k > nl by testing (13 . 18[) 
by (ji — k) + . In the present case, the equation to be tested is (I3.24p instead of (13.181) . 
and a more elaborate procedure is needed. First of all, we check that (fi — k) + is an 
admissible test function (this is not obvious since V/i might not exist in the usual sense). 
We recall that A is a strictly increasing mapping from [0, +00) onto itself and that A -1 is 
Lipschitz continuous on the interval [s*, +00), where s* := A'(r*), due to (13.41) . Therefore, 
we can choose a strictly increasing map A* : [0, +00) — > [0, +00) that is globally Lipschitz 
continuous and coincides with A -1 on [s*, +00). Hence, we have A„(A(r)) = r for every 
r > r* and A*(A(r)) < for r < 7%. It follows that (r — k) + = (A*(A(r)) — k) + 
for every r > if k > r*. On the other hand, K t (K(ji)) G H l {fd,T;H) n L 2 (0,T;\/) 
by (13.221) . Hence, (fi — k) + enjoys the same regularity and is an admissible test function 
in (I3.24p provided that k > r*. Thus, we assume k > max{/iQ, r*} from now on. We have 
from fl3T24"j) 

[\d t [(1 + 2<7(p))//] (s), (ji(s) - k) + ) ds+ [ VKQi) ■ V(ji - k) + 

JO JQt 

^dtg(p) (fJ. - k) + 
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for every t £ [0, T], and a simple rearrangement yields 

\d t [(l + 2g(p))(p-k)](s),(p(s)-k) + )ds + [ VK(fi) • V(/i - k) + 

JQt 

= [ d t g(p)\(p-k) + \ 2 -k [ d t g(p) 0* - k) + . (5.3) 

JQt JQt 

Now, noting that l/(l+2#(p)) £ H\0,T; V)nL°°(0,T; W) by (l3~T2l and our assumptions 
on g (recall (I3.6l) -( l3~8l) ). we apply Lemma ISTTl with w := (1 + 2g(p))(p — k) and 7 := 
1/(1 + 2g(p)) and transform the first term on the left-hand side as follows: 



(d t [(l + 2g(p))(p-k)](s),(p(s)-k)+}ds= / (p-k)d t [(l + 2g(p))(p-k) + ] 

JQt 

2d t g{p)\(p-k) + \ 2 + [ (p-k)(l + 2g(p))d t (p-k) + 



JQt 

d t [(l + 2g(p))\(p-k) + \ 2 ] + I d t g{p)\(p-k) + 



1 

~ 2 

On the other hand, we have a.e. in the set where p > k 



v(/i - k) + = v/i = v^^w) = (K- 1 )'(K( / i))VK( yU ) = -jL vif a*). 

K(/i) 

Finally, (yu(0) — /c) + = since k > p$. Therefore, ( 15 .31) becomes 

\ hl + 2g(p(t)))mt)-k) + \ 2 + [ K(p)Mp-k) + \ 2 = -k [ d t g(p) (p - fc)+ . 

.In ./r>. ./r>, 



Since g is nonnegative and n(r) > for r > k (because k > k*), we obtain 
\(p(t) - k) + \ 2 + kJ \V(p-k) + \ 2 <k I \d t g{p)\{p-k) + . 



1 

2 

At this point, the argument used in j8] can be repeated without changes, essentially 
Indeed, the analog of (13. 14ft is never used there, and the whole proof is based just on the 
regularity d t p £ L°°(0,T; H) n L 2 (0,T; V). In the present case, we have to exploit the 
same regularity for dtg(p) which follows from (I3.12p . □ 

Remark 5.2. As observed in Remark 13.41 the component p of any weak solution is a 
weakly continuous if -valued function. We show that 

p £ L7°([0,T];L p (fi)) for p£ [1,2), and for p £ [1, +00) if p Q £ L°°(n). (5.4) 

Assume t n — > t £ [0,T]. We prove that p(t n ) — > p{t) strongly in L P (Q) for p like 
in (15 .4p . In fact, p(t n ) is bounded in L 2 (Q) in the general case, and in L°°(Q) if /io is 
bounded (thanks to Theorem 13. 7p . So, the desired convergence is proved once we show 
that p{t n ) — > p(t) a.e. in Q, at least for a subsequence. We observe that (I3.22p implies 
K(p) £ C°([0, T]; if ). Hence, K(p(t n )) — > K(p) a.e. in fl, at least for a subsequence. As 
i\~ _1 is continuous, the claim follows. 
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Proof of Theorem \3.8l As ( I3.26P holds, g(p) is continuous and g'(p) is bounded. On the 
other hand, p is bounded too (by Theorem 13.71 since p is bounded). Hence, (13 . 1 8[) can 
be seen as a linear uniformly parabolic equation for p with continuous coefficients and a 
right-hand side belonging to L°°(0,T; H). We have indeed 



By owing to p G W and optimal L p -L 9 -regularity results (see, e.g., [131 Thm. 2.3]), we 



Remark 5.3. We notice that the same result (13.271) holds under an assumption on p 
that is weaker than po G W. The optimal condition involves a proper Besov space and 
can give a similar result for a fixed p. We are going to use (13.271) just with p = 4 in our 
proof of the uniqueness of the solution. It follows that uniqueness still holds for a less 
regular p . 

Remark 5.4. We observe that the case corresponding to an empty interior of D((3) is 
completely trivial. Indeed, if D(j3) = {ro}, then p takes the constant value ro, p solves the 
corresponding heat equation, and £ is computed from (13.141) . In the opposite case, further 
regularity can be proved under suitable assumptions on the initial data. For instance, by 
supposing p to be bounded and nonnegative, we note that (13.141) yields 



So, by assuming that infp an d sup po belong to the interior of D((3), one can easily 
derive that £ G L°°(Q). Indeed, one can formally multiply by |£| p-1 sign£ and estimate 
|| £ ||p uniformly with respect to p if this assumption on po is satisfied. This implies that 
p G H/ 1 ' p (0,T;L p (fi)) n U>(0,T;W 2 ' p (tt)) for every p < +oo whenever p is smooth 
enough. However, no further regularity can be proved, in general, since (13.141) cannot 
be differentiated, unless f3 is particular, e.g., like in [8]. By the way, in that case, the 
condition £ G L°°(Q) is equivalent to inf p > and sup p < 1. More generally, if D(f3) 
is an open interval (a, b) and (3 is a smooth function, then £ = /3(p) G L°°(Q) implies 
that inf p > a and sup p < b and a bootstrap technique using both equations can lead to 
higher regularity. 

6 Uniqueness 

In this section, we prove Theorem 13.91 We observe that the uniqueness of the third 
component £ follows by comparison in (13.141) once we prove that the pair (p, p) is unique. 
So, we deal with the first two components, only, and remind the reader that we can use 
the further regularity given by Theorem 13.81 In particular, by accounting also for (13 . 1 2[) 
and (13.281) . we have 



d t p - (1 + 2g(p)) 1 Ap = - (1 + 2g(p)) ~ l p g'{p) d t p. 



infer that fl3T27]) holds. 



□ 



d t p - Ap + £ = pg'(p) - 7r(p) G L°°(Q). 



\Vp\ G L 4 (0 



T;L 6 (fi)) and | Vp| G L 4 (0, T; L 6 



(«)) 



(6-1) 



for every solution. First of all, we rewrite equation (13.181) in the form 



d t (p/a(p)) - a(p)Ap = 0, 



(6.2) 



21 



Existence - uniqueness for a Cahn-Hilliard system with viscosity 



where the function a : [0, +00) — > (0, +00) is defined by 

a{r) := (l + 2g(r))~ 1/2 for r > 0. (6.3) 

More precisely, we consider the variational formulation of (16. 2p that accounts for the 
homogeneous Neumann boundary condition and involves a related unknown function, 
namely 



z: " 



and I d t z(t)v+ [ V(a{p{t)) z(t)) ■ V(a(p(t))v) = 
Jn Jn 



a(p) 

for a. a. t G (0,T) and for every v G V. (6-4) 

Notice that z is bounded since both p and p are. Indeed, (I3.26P holds and Theorem 13.71 
can be applied since W C L°°(Q). Moreover, z satisfies the analogue of (16. ip . At this 
point, we pick two solutions (yUj,pj,£j), i — 1,2, and set 

:= a(pi) and Zi := Pi/ai for i = 1,2 

so that (zi, Pi) satisfy (16 .4p . In the subsequent estimates the (varying) value of the constant 
c may even depend on the considered solutions, e.g., through ||zi||oo- Our method proceeds 
as follows. We write (16 .4p for both solutions and choose v = z\ — z<i in the difference. Then 
we integrate over (0,t), where t e (0, T) is arbitrary. At the same time, we write (I3.14p 
for both solutions and multiply the difference by p± — p2- Then we integrate over Q t . 
Finally, we take a suitable linear combination of the resulting equalities and perform a 
number of estimates that lead us to apply the Gronwall lemma. However, in order to 
simplify notation and make the proof more readable, we set 

p:= pi- p 2 , p:= pi- P2, f := £1 - 6, z\=z\ — z 2 , and a := a x - a 2 , 

and prepare some auxiliary material before starting. The next inequalities account for 
the boundedness and the Lipschitz continuity of a, a', and 1/a on the range of p (recall 
that p is bounded). We have 

a\ = \a(pi) - a(p 2 )\ < c\p\, 

Va\ = \a\ Pl )Vp + (a'(pi) - a'(p 2 ))Vp 2 \ < c\Vp\ + c\ Vp 2 | \p\, 
Va" 1 ) < c|Vpi|, 

p\ < I aj \z 2 \ + a 2 \z\ < c\a\ + c\z\ < c\p\ + c\z\, 
\Vz\ = \V(a^ 1 (a 1 z))\ < c|V(ai^)| + c\Vpi\ \z\. 

In what follows, we will repeatedly use these inequalities without reminding the reader. 
Lemma 6.1. We have, for every t G [0,T], 

/ \Wz\ 2 <c[ |V(M| 2 + c f\l + \\Vp 1 (s)\\i)\\z(s)\\ 2 2 ds. (6.5) 
JQt JQt Jo 

Proof. By the preliminary inequalities just stated, we have 

\Vz\ 2 <cf \V(a lZ )\ 2 + c[ |Vpi| 2 |z| 2 . (6.6) 
Qt JQt JQt 
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Now, by using the Holder and Sobolev inequalities (see (I3.28P ). we obtain 

|Vpi| 2 M 2 < fllVp^wlMs^Wz^hds 







<c ||V Pl ( S )|| 2 (||V,( S )|| 2 + P( S )|| 2 )||,( S )|| 2 ^ 
Jo 

<c /W(*)ll«ll*(*)ll2<fo 

Jo 

+ e [ \Vz\ 2 + c e Avpi(s)ll3ll*(s)ll> 

JQt Jo 

<e [ \Vz\ 2 + c £ [\l + \\V Pl (s) 
JQt Jo 



18 \W)\\lds, 



where e > is arbitrary. Hence, (16.51) follows by combining this with (16.61) and then 
choosing e small enough. □ 

Lemma 6.2. Let k e L 4 (0, T; L 6 (fi)). Then we have 

I k 2 (\z\ 2 + \p\ 2 )<s j (|V( ai z)| a + |Vp| a ) 

JQt JQt 

+ c £ f\l + \\V Pl (s)\\t + \\k(s)\\t){\\z(s)\\l + \\p(s)\\l) ds (6.7) 
Jo 

for every e > and every t G [0, T}. 

Proof. By the Holder and Sobolev inequalities (see (I3.28|) ). we have 

k 2 (\z\ 2 + |p| 2 ) < fwmWl (\\z(s)\\ 6 \\z(s)\\ 2 + ||p( a )|| 6 \\p(s)\\ 2 ) ds 
Jo 

< c /W)ll6(l|V^)|| 2 \\z(s)\\ 2 + ||,( S )|| 2 + ||VpOO|| 2 ||p( S )|| 2 + \\p(s)\\l) ds 
Jo 

< e [\\\Vz(s)\\l + \\Vp(s)\\l) ds + c £ [\\\k(s)\\t + \\k(s)\\t) {\\z(s)\\j + \\p(s)\\l) ds. 



By applying Lemma 16.11 and denoting the constant that appears in ( 16. 5ft by C, we can 
continue and obtain 



k 2 (\z\ 2 + \ P \ 2 ) 

<s(c f \V(a lZ )\ 2 + C [\l + \\V Pl (s)\\t)\\z(s)\\lds+ [ \V P \ 2 ) 

\ JQt Jo JQ t J 

+ c £ [\l + \\k(s)\\t)(\\z(s)\\l + \\p(s)\\l)ds. 
Jo 

Hence, (16. Tp immediately follows. □ 
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At this point, we can start with our program. However, in order to make the argument 
more transparent, we deal with the first equation only, for a while. We have 

/ W)\ 2 + I (V(aizi) • V(auz) - V{a 2 z 2 ) • V{a 2 z)) = 0. 

1 Jn JQ t 

It is convenient to transform the last integrand as follows: 
V(aiZi) ■ V(aiz) - V(a 2 z 2 ) ■ V(a 2 z) 

= |V(ai2)| 2 + V(aiz 2 ) • V(ai^) - V(a 2 z 2 ) • V(ai*) + V(a 2 z 2 ) • V(az) 
= |V(ai2)| 2 + V(az 2 ) • V(aiz) + V/i 2 • V(az). 

Then, the above equality becomes 

\z{t)\ 2 + I \V(a lZ )\ 2 = - [ V{az 2 )-V{a x z)- ! Vp 2 -V(az), (6.8) 

JQt JQt JQt 

and we estimate each term of the right-hand side separately. We immediately have 

-/ V(az 2 ) ■ V(a lZ ) < - [ \\7(a lZ )\ 2 + 2 [ (z 2 2 \ Va| 2 + a 2 \ Vz 2 | 2 ) 
JQt 4 JQt JQt 

<\f \V(a lZ )\ 2 + c[ {\Vp\ 2 +\Vp 2 \ 2 \p\ 2 )+c [ \Vz 2 \ 2 \p\ 2 

4 JQt JQt JQt 

<- f Mayztf + df |Vp| 2 + c/ (|Vp 2 | 2 + |Vz 2 | 2 )|p| 2 , (6.9) 

4 JQt JQt JQt 

where we have marked the constant that we want to refer to by terming it C\. We treat 
the last term of ( 16. 81) as follows: 

-/ VpL 2 -V(az)< |VAt 2 |(|a| |Vz| + \z\ |Va|) 

JQt JQt 

< c f |V^|(|V( ai z)| |p| + \z\ \V Pl \ \p\ + \z\ |Vp| + \z\ |Vp 2 | |p|) 

JQt 

<- f |V(M| 2 + C / |Vp 2 | 2 |p| 2 + / |v P | 2 

4 JQt JQt JQt 

+ cf \Vp 2 \ 2 \z\ 2 + c f (|Vp!| 2 +|Vp 2 | 2 )|p| 2 . (6.10) 

JQt JQt 

Now, we deal with the second equation. Testing the difference of ( I3.14p by p as mentioned 
at the beginning, easily yields 

\ [ [ |Vp| 2 + ( Zp= ( (pi^(pi)-p 2 ^(p 2 )-7r(pi) + 7r(p 2 ))p. (6.11) 

A Jn JQ t JQt JQt 

We note that the last integral on the left-hand side of (16. lip is nonnegative by mono- 
tonicity, while the integrand on the right-hand side can be estimated as follows: 

(p\g'{Pi) - lMt$f(p%) ~ 7r(pi) + 7r(p 2 ))p 

< (H W(pi)\ + N \g\ P i) - g'( P 2)\ + Hpx) - n{ P2 )\)\p\ 

< \g\Pi)\ \P\ \P\ + c\p 2 \ |p| 2 + c|p| 2 < c(|/i| 2 + |p| 2 ) < c(\z\ 2 + |p| 2 ). 
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Thus, by inspecting the coefficients of the integral jg t | Vp| 2 that appear on the right-hand 
sides of (I6.9p and (I6.10p . it is clear that it is convenient to multiply (16. lip by G\ + 2 
before adding it to (16. 8p . Once such a care is taken, it is straightforward to deduce that 

/ \ z (t)\* + I |V(M| 2 + / | P (t)| 2 + / |v P | 2 

Jn JQt Jn JQ t 

<c[ (|V/i 2 | 2 + l) \z\ 2 + c [ (|V/i 2 | 2 + |Vpi| 2 + |Vp 2 | 2 + |V2 2 | 2 + 1) |p| 2 

JQt JQt 

<c (|V/i 2 | + |Vpi| + |Vp 2 | + \Vz 2 \ + l) 2 (\z\ 2 + |p| 2 ). 

JQt 

At this point, we recall that (16.11) hold for p i; p i; and Zj, so that we can apply Lemma [6721 
with k = |Vp 2 | + |Vpi| + |Vp 2 | + |Vz 2 | + 1. Then, we choose e > small enough and 
use the Gronwall lemma. We conclude that z = and p = 0, whence (pi,pi) = (p 2 ,p 2 ) 
follows. 
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